Interaction effects on the conductance peak height statistics in quantum dots 
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A random interaction matrix model is used to study the 
statistics of conductance peak heights in Coulomb blockade 
quantum dots. When the single-particle dynamics conserves 
time-reversal symmetry, the peak height statistics is insensi- 
tive to the interaction strength. But when the single-particle 
dynamics breaks time-reversal symmetry, the peak height 
statistics exhibits a crossover from unitary to orthogonal sym- 
metry as the interaction strength increases. This crossover is 
driven by the time-reversal symmetry of the interaction. Our 
random interaction matrix model describes features of both 
the measured peak height and peak spacing statistics. 

PACS numbers: 73.23.Hk, 05.45+b, 73.20.Dx, 73.23.-b 



Random matrix theory (RMT) provides a useful tool 
for describing the universal statistical fluctuations of the 
spectrum and eigenfunctions of a quantum system whose 
associated classical dynamics is chaotic. RMT has been 
successfully applied to the study of mesoscopic phenom- 
ena in quantum dots - submicron 2D devices where elec- 
trons are confined by electrostatic potentials In 
dots with irregular shapes the single-electron dynamics 
is mostly chaotic, and the use of RMT is justified within 
a single-particle framework. In dots that are strongly 
coupled to leads, or open dots, the approximation of non- 
interacting quasi-particles is reasonable, and interactions 
can be considered indirectly through their effect on the 
electron coherence time. When a finite dephasing time is 
included, RMT can describe quantitatively the observed 
conductance statistics in open dots . 

As the dot-leads coupling is made weaker, the charge 
on the dot becomes quantized, and electron-electron in- 
teractions cannot be ignored. In such almost closed dots, 
the conductance displays sharp peaks versus the gate 
voltage. The measured peak height distributions j|]4[] 
were found to agree with RMT predictions M. However, 
it is not clear why RMT should describe the peak height 
statistics in dots with strong electron-electron interac- 
tions. For example, the spacings between peaks were 
observed to have a Gaussian-like distribution and 
not the Wigner-Dyson distribution that is expected in 
the constant interaction model plus single-particle RMT. 
Numerical simulations of an Anderson model of a small 
disordered dot (< 10 electrons) with Coulomb interac- 
tions also found Gaussian peak spacing distributions || , 
while the peak height distributions showed only a weak 
dependence on interactions M. 



RMT is not limited to non-interacting systems. In- 
deed it was first introduced to explain the neutron res- 
onance data in the compound nucleus. But the neutron 
resonances are measured at finite excitations, while the 
linear conductance experiments in quantum dots probe 
the ground state of the system with a varying number 
of electrons. To understand the statistics that emerges 
from the interplay between one-body chaos and interac- 
tions in these dots, it is necessary to construct a ran- 
dom matrix model that includes interactions explicitly 
and reduces to one-body RMT in the absence of inter- 
actions [|J . We can break the time-reversal symmetry of 
the one-body Hamiltonian (e.g., with a magnetic field), 
but the time-reversal symmetry of the two-body inter- 
action should be preserved. An important question is 
whether such a random matrix model can reproduce both 
the measured peak height and peak spacing statistics. 

A random interaction matrix model (RIMM) was in- 
troduced in Ref. [jl0| to study the peak spacing distribu- 
tion in chaotic dots. The model combines a random two- 
body interaction pd[ | with a random one-body Hamil- 
tonian, and describes a crossover of the peak spacing 
statistics from a Wigner-Dyson distribution to a Gaus- 
sian distribution in terms of a parameter that measures 
the fluctuations of the interaction matrix elements. Here, 
we use the RIMM to study the effects of interactions on 
the conductance peak height statistics. For a GOE one- 
body statistics, we find that the partial width (to decay 
into one of the leads) has a GOE Porter-Thomas distri- 
bution independently of the interaction strength. How- 
ever, for a GUE one-body statistics, we find a crossover 
from a GUE to a GOE Porter- Thomas distribution as 
the interaction strength increases. It is well known that, 
in the absence of interactions, a complete crossover from 
GOE to GUE statistics can be induced by a magnetic 
field. Our results indicate that this crossover is not com- 
plete when strong interactions are turned on because of 
a competition between an asymptotic symmetry of the 
one-body Hamiltonian and the GOE symmetry of the 
time-reversal-conserving two-body interaction. Finally, 
for the case of GUE one-body statistics, we find that a 
Gaussian-like peak spacing distribution develops at in- 
teraction strengths for which the peak height statistics is 
still close to GUE. This explains the observed peak spac- 
ing statistics |||7| and peak height statistics within 
a single random matrix model. 

The n-th conductance peak height G n corresponds to 
the tunneling of an electron into a dot with n— 1 electrons 
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to form a dot with n electrons. 
G n = (e 2 //i)(7rf/4fcT) 5 „, where 
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9n = 



n n 



At low temperatures, 
(1) 



T 1 ^ is the partial width of the ground state of the n- 
electron dot to decay into an electron in the left (right) 
lead and the ground state of the dot with n — 1 electrons: 



($ g . s .(n)| ^(r)\ <D g 



(2) 



ijv (r) is the creation operator of an electron at the point r 
[r = Tim for the left (right) point contact], and <& gs .(rt) 
is the ground state wavefunction of the n-electron dot. 
For non-interacting electrons, $ g . s .(n) is a Slater deter- 
minant of the n lowest single-particle eigenfunctions in 
the dot, and Eq. (||) reduces to T n oc \(f> n (r)\ 2 , where <j) n 
is the n-th single-particle wavefunction. If the single- 
particle dynamics is chaotic, |</>„(r)| 2 satisfies Porter- 
Thomas statistics, leading to universal distributions of 
the conductance peak heights (|l|) that are sensitive only 
to the underlying symmetry class ||. 

To determine how interactions might modify the 
Porter-Thomas statistics of the partial widths, we use 
the RIMM of spinless interacting fermions fiof] 
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The one-body matrix elements hij are chosen from the 
appropriate Gaussian random matrix ensemble, while 
the anti-symmctrized two-body matrix elements Uij-ki 



Uij-kl 



:-ik form a GOE in the two-particle space |TT| 



P(h) 



P(u)c<e- Trfi2 / 2C/2 



(4) 



The states \i) = a]|0) describe a fixed basis of m single- 
particle states, h is an m x m GOE (GUE) matrix when 
the single-particle dynamics conserves (breaks) time- 
reversal symmetry. The two-body interaction is assumed 
to preserve time-reversal symmetry and forms a GOE, ir- 
respective of the symmetry of the one-body Hamiltonian. 

To calculate the statistics of the partial widths T„, 
we expand tp\r) — J2ii J i( r ) a li where -0i( r ) = (r\i) is 
the wavefunction of the fixed state It follows from 
the orthogonal invariance of the ensemble (^) that the 
statistics of T n is identical to the statistics of 



(<£ g . s .(n)| a] |$ g . s .(n-l)) 



(5) 



for any i. 

For each realization H of the ensemble (pj), we calculate 
the ground states for n — 1 and n electrons, and compute 
rjj using (^) . The distributions of the normalized width 
r = r/r for a GOE single-particle statistics are shown in 



Fig. |l| for several values of U/ A. We show P(ln f ) rather 
than P(r) in order to display more clearly the small val- 
ues of f. The distributions are independent of U/A and 
are well described by the GOE Porter-Thomas distribu- 
tion P GOE (lnf) = (f/2) 1 / 2 exp(-f/2) (solid line). As a 
reference we also show the GUE Porter-Thomas distribu- 
tion PcuE(lnr) = Texp(— f) (dashed line). The conduc- 
tance peak heights are calculated from (Q) using two un- 
corrected "sites" i and j for the left and right leads. The 
peak height distributions P(g) (for a GOE h) are shown 
in the insets of Fig. [l]. They are all in good agreement 
with the GOE distribution P G OE(.g) = ^2/nge- 2g (solid 
lines) irrespective of U/A. The dashed lines describe the 
GUE distribution P(g) = 4g[K (2g) + J ftT 1 (2 5 )] e - 29 (K 
and K\ are Bessel functions). 
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FIG. 1. Width and conductance peak height distributions 
in the RIMM with GOE one-body statistics. P(ln Y) us shown 
versus lnf (f is the normalized width (&)) for m = 12, n = 4 
and U/A — (circles), 2.4 (squares), 4 (diamonds) and 8 
(triangles). The solid and dashed lines are the GOE and 
GUE Porter- Thomas distributions, respectively. Insets: the 
peak height distributions P(g) versus g in a log-log scale (left 
inset) and in a linear scale (right inset) for the same cases 
shown in the main figure. The solid and dashed lines are the 
GOE and GUE peak height distributions, respectively. 

The results for a GUE one-body statistics are shown in 
Fig. ^ for the same values of U/A as in Fig. When the 
interaction strength increases, the width and peak height 
distributions make a crossover from the corresponding 
GUE distribution (at U = 0) to the GOE distribution (at 
large U). Equivalcntly, the transition from GOE to GUE 
statistics due to a time-reversal symmetry-breaking one- 
body field is not complete because of the competing GOE 
symmetry of the two-body interaction. The crossover dis- 
tributions are compared with distributions obtained from 
an RMT ensemble that describes the crossover between 
the orthogonal and unitary symmetries. This ensemble 
is H = S + iaA , where S and A are N x N symmetric 
and antisymmetric uncorrelated Gaussian matrices and 
a is a parameter jl3| . Its wavefunction statistics depends 
on the parameter A = a^/N /n. In particular, the width 
distribution is given in closed form by [pT|Jl^] 
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p x (t) = j dtp x (ty / 



1-t 4 



(6) 



where Jo is the modified Bessel function of order zero. 
t in (H) is a "shape" parameter that fluctuates in the 
interval [0, 1] according to a known distribution |l4|-|l(| 

Px(t) = tt 2 (1A 3 - t) A^-^^Ct-V*) 2 {M x)+ 

^(i + ^ 1 ) 2 /4-l/(2 7 r 2 A 2 )] [l-^(A)]} , (7) 

where <pi(X) = J* dy e -2 Tr 2 A 2 (i- y 2 )_ Eq ^ re duces to 
the GOE and GUE Porter-Thomas distributions for A = 
and A> 1, respectively. In practice, the crossover from 
the GOE to the GUE already occurs for A ~ 1. 




FIG. 2. Width and conductance peak height distributions 
in the RIMM with GUE one-body statistics. P(ln f) is shown 
versus In f for the same cases as in Fig. []j The short-dashed 
lines are fits to the theoretical distribution (^) describing the 
RMT crossover from GUE (A = oo, dashed line) to GOE 
(A = 0, solid line) with A = 0.28, 0.17 and 0.08, respectively. 
Insets: P(g) in a log-log scale (left inset) and in a linear scale 
(right inset) for the same cases shown in the main figure. The 
short-dashed lines are the analytic peak height distributions 
in the crossover regime |l^| between the GUE (dashed) and 
GOE (solid) distributions. 

For each U/A, we find A by fitting the distribution 
(^) to the computed width distribution. The distribu- 
tions ([}]), shown by the short-dashed lines in Fig. g, ac- 
curately describe the width distributions of the RIMM 
with a GUE h. Good agreement with closed RMT ex- 
pressions is also obtained for the peak height distri- 
butions P\(g) shown in the insets of Fig. g. 

An important issue is the universality associated with 
the RIMM (§) and ©■ The model depends on three 
parameters: m, n, and U/A. The top panel of Fig. || 
shows, for a GUE one-body h, the crossover parameter 
A as a function of U/A for m = 10 and n = 4,5,6,7 
(symbols), and for a "reference" case m = 12 and n = 4 
(solid line). The curves depend on both m and n, but 
can all be scaled on the reference curve after scaling the 
interaction strength by a constant, U e g = f(m, n)U/A 



(see bottom panel of Fig. [|). The values of the scaling 
factor /(m, n), shown in the bottom inset of Fig. ||, are 
essentially the same as those needed to obtain a universal 
peak spacing statistics pToj ] . We conclude that the peak 
spacing statistics as well as the peak height statistics in 
the RIMM depend only on a single parameter U e g. 
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FIG. 3. Top: The RMT crossover parameter A as a func- 
tion of U/A for m = 10 and n — 4 (circles), 5 (squares), 6 (dia- 
monds) , and 7 (triangles) . The solid line is the reference curve 
m = 12, n = 4 Top inset: The average width f (U)/f(0) ver- 
sus U/A for the same cases shown in the main figure. The 
lines are fit to f (U) - f (00) = (f (0) - f (oo))/(l + b U 2 /A 2 ). 
Bottom: The curves from the top panel scale on the reference 
curve (solid line) when plotted versus U e s = f(m,n)U/A. 
The bottom inset shows the scaling factor f(m, n) (open sym- 
bols) versus n for m = 10 (circles), 12 (squares), and 14 (dia- 
monds). The scaling factors are nearly the same as the ones 
derived from the peak spacing statistics (solid symbols). 

In the range U c g < 1, A > 1 and the statistics is essen- 
tially GUE. For Ucs ~ 1 — 1.5, the peak height statistics 
is close to GUE while the peak spacings already follow a 
Gaussian distribution. This explains the measured RMT- 
like peak height distributions and the Gaussian-like 
shape of the peak spacing distributions ||,|7j within a sin- 
gle random matrix model. We remark that the small de- 
viations from GUE statistics observed in the experiment 
of Ref. H and in the calculations of Ref. || in the pres- 
ence of a magentic field are consistent with a crossover 
from GUE to GOE. 

The average width T is a monotonically decreasing 
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function of U/A and saturates at large values of U. The 
top inset of Fig. || shows T(U) /f (0) as a function of U/ A 
for several m and n. This dependence is well described 
by f (17) - f(oo) = (f (0) - f(oo))/(l + bU 2 /A 2 ). The 
parameter r(oo) depends on m and n, while & is found 
to be independent of m and to depend only weakly on n. 

Next we compare the predictions of the RIMM with 
those of a model of a quantum dot. We studied a 2D 
Anderson model with on-site disorder parameter W and 
hopping matrix element V = 1. The electrons are inter- 
acting with a Coulomb interaction whose strength over 
one lattice spacing a is U c = e 2 /a B. We choose peri- 
odic boundary conditions in both directions so that the 
average width and width statistics are independent of 
the specific sites to which the leads are attached. With 
these boundary conditions, a background charge (that 
was included in the calculations of Rcf. ||) only shifts 
the energy levels by an rt-dependent constant but does 
not affect the wavef unctions. 




FIG. 4. Width distributions in a 4 x 5 Anderson model 
plus Coulomb interactions. Left: P(lnP) in the absence of 
a magnetic field and for several values of the Coulomb inter- 
action strength: U c = (circles), 2 (squares), 4 (diamonds) 
and 6 (triangles). The disorder parameter is W = 5. The 
solid (dashed) line is the GOE (CUE) Porter-Thomas distri- 
bution. Right: P(lnP) for U c — (circles) and 12 (triangles) 
in the presence of time-reversal symmetry-breaking magnetic 
flux $ = 0.14$o and for a disorder strength of W = 3. The 
short-dashed line is the distribution (^) with A = 0.17. 

The left panel of Fig. |J shows the width statistics in 
the absence of magnetic flux for a 4 x 5 lattice with 
disorder parameter W — 5 and for n — 4 electrons. 
The distributions are approximately described by the 
GOE Porter-Thomas distribution (solid line) irrespective 
of the value of U c and in agreement with the RIMM. 
The dependence of the average width T on the interac- 
tion strength (not shown) is similar to that observed in 
the RIMM. The right panel shows width distributions 
in the presence of magnetic flux <E> = 0.14$o and for 
W = 3. The U c = distribution (circles) agrees with 
the GUE Porter-Thomas distribution (dashed line), while 
the U c — 12 distribution (triangles) is described by (0) 



with A = 0.17. This is the distribution obtained in the 
RIMM for U c s m 4. For weaker interactions the calcu- 
lated width distributions exhibit some deviations from 
(^). We note however that, for the small lattices used, it 
is difficult to find a disorder strength for which the model 
is in the metallic diffusive regime and displays universal 
RMT statistics. In particular, in the presence of a mag- 
netic flux we could not find values of W for which both 
the spectral and wavefunction statistics are GUE. 

In conclusion, we have investigated the width and peak 
height statistics in Coulomb-blockade quantum dots us- 
ing a random interaction matrix model with interactions 
that preserve time-reversal symmetry. For a GOE one- 
body symmetry the statistics is insensitive to the inter- 
action strength. However, at strong interactions, a time- 
reversal symmetry-breaking field leads only to a partial 
crossover from GOE to GUE statistics. Our random in- 
teraction matrix model can reproduce both the observed 
Gaussian-like shape of the peak spacing distribution and 
the RMT statistics of the peak height distributions. 
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